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Abstract 

We study a spin system on a large box with both Ising interaction 
and Sherrington-Kirpatrick couphngs, in the presence of an external 
field. Our results are: (i) existence of the pressure in the limit of an in- 
finite box. When both Ising and Sherrington-Kirpatrick temperatures 
are high enough, we prove that: (ii) the value of the pressure is given 
by a suitable replica symmetric solution, and (iii) the fluctuations of 
the pressure are of order of the inverse of the square of the volume 
with a normal distribution in the limit. In this regime, the pressure 
can be expressed in terms of random field Ising models. 
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1 Introduction 

We consider a (i-dimensional Ising model in a magnetic field h, perturbed by 
a mean field interaction of spin-glass type. The Hamiltonian contains two 
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parameters, (3 and k, which play the role of two inverse temperatures. When 
P = the model reduces to the Ising model at temperature while for k = 
one recovers the Sherrington-Kirkpatrick (SK) model at temperature 
The understanding of the SK model has recently witnessed great progress 
(see, e.g., [Hj ^M)- The main interest in the analysis of this model is 
the possibility of investigating the robustness of the phenomena typical of 
mean field spin glass models, in the presence of additional interactions of 
non-mean-field character. 

The model has been previously considered in [5] under the additional 
assumptions that h = and that the Ising model is ferromagnetic. Under 
these conditions it was proven that, if (3 and k are small enough, the infinite 
volume pressure is given by the sum of the Ising pressure and of the SK one, 
at the respective temperatures. Moreover, the disorder fluctuations of the 
pressure were found to be of order 1/V, V being the volume of the system, 
and to satisfy a central limit theorem. 

In the present paper, the Ising interaction decays exponentially fast with 
distance, but is not necessarily ferromagnetic. It turns out that the presence 
of the magnetic field changes quahtatively the picture with respect to jS]. 
Indeed we find, still for k and (3 small enough, that the limit pressure is given 
in terms of the pressure of an Ising model with random external field, the 
strength of the randomness being related to the typical value of the overlap 
between two replicas of the system. It is remarkable, though natural, that 
the random field Ising model, which has its own interest [2 , plays such an 
important role in our model. Note also that the pressure of our model can be 
computed only via the thermodynamic limit of another disordered system. 
This is contrast with the case h = we mentioned above, and also of course 
with the case n = of the standard SK model. The second difference is 
that the fluctuations of the pressure in presence of h satisfy a central limit 
theorem on the scale 1/VV rather than 1/V. The same phenomenon is 
known to happen in the case of the usual SK model, see for instance jTj |3] 

In the region of thermodynamic parameters we consider, the system is 
in a "replica-symmetric" (RS) phase, the overlap between two independent 
replicas being a non-random value in the thermodynamic limit. Our methods 
fail beyond some values (3o{h) and Ko{h), which we believe to be an artifact of 
our approach, rather than representing the true boundary of the RS region. 
The same inconvenient has been previously encountered in the analysis of 
the SK model ^21; [HI- In principle one could improve these values by em- 
ploying the replica symmetry breaking scheme which for instance enabled 
M. Talagrand fH] to control the whole replica-symmetric region of the SK 
model, and later the entire phase space ^1] . However, in the present case (as 
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well as in |3j) one of the reasons why we do not reach the true critical line 
is due to an incomplete control of the underlying random field Ising model, 
and this problem would not be fixed by the methods of |16^ . For this reason, 
we prefer to use a generalization of the technically simpler "quadratic replica 
coupling" technique introduced in pi. 

It would be of course an interesting challenge to go beyond the present 
approach, and to deal with lower-temperature situations where the Ising-SK 
system possibly shows a RSB-like behavior. 



2 Description of the model and results 

The model we consider is defined on the d-dimensional hypercubic box An = 
{—N, ■ ■ ■ , A^}'^ and its partition function is: 

Z^in, P, h-J)= Yl ( - PH^^'i^-^ J) + hJ2^i]- (1' 

o-e{_i,+i}Ajv \ ieAjv / 

The Hamiltonian of the Sherrington-Kirkpatrick (SK) model is defined as 



2lA7v|.tr '''''''' 



and the couplings Jjj are i.i.d. Gaussian random variables A/'(0, 1). On the 
other hand, the Hamiltonian of the Ising model is 



2 

«j'eAjv 

where we assume that the interaction decays exponentially, i.e., 

\K{{}\<C,e-^'\'\ (2) 

for some Ci,C2 > 0. We do not require the interaction to be ferromagnetic. 
The finite volume (disorder-dependent) pressure is defined as usual as 

PN{K,P,h; J) = -r^—r log ZN{K,P,h; J). 

Later, we will need to consider the (Gaussian) Random Field Ising Model, 
defined by the partition function 



Zn {i^,h,'J,J) 



J2 exp I -KHl,{a) + E a,{h + 7 J,) j 
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Jj being i.i.d. standard Gaussian variables, and being also independent from 
the Jjj's in the sequel. The existence of the infinite volume pressure of the 
RFIM, 



p™(«:,/i,7) = lim EpT''i'^,h,r,J) 

,RFIM 



= a.s. — lim {K,h,^; J) , 

N^oo 

is a well known consequence of additivity and of the ergodic theorem, e.g. 

Our main results can be summarized as follows: 
Theorem 1. For all /i G M, and all K,,f3 > 0, the limit 



exists and 



p{K,P,h)= lim Kpn{k, P, h; J) (3) 

TV— >oo 



PN{K,f3,h;J)^-^p{K,(3,h) (4) 



for a.e. J and in the U'-norm (p G [1, oo) ). 

Theorem 2. For all /i G M, there exist Ko{h) > and (3o{h) > such that, 
forO< K< Ko{h) andO <P < po{h) 

p{n,P,h)= ini (p^^^^(n,h,P^) + ^{l-qA (5) 
o<i3<i \ 4 y 

In the = 0, it is established in [7j that the infimum in (jS)) is 

achieved at a unique point. As explained below (see Section IT^ . we take 
small enough kq^H) so that, for < k < Kq^H) and < (3 < /3o{h), the 
infimum in is achieved at a unique q. In the following, we will always 
denote by g = f3, h) the value that realizes the infimum in 

We emphasize that no{h) is taken small enough so that the RFIM is 
inside Dobrushin's uniqueness region for every realization of the external 
fields Ji,i G Ij^. Let be the unique infinite volume Gibbs measure 

for the RFIM with 7 = j3y/q, depending on Jj, i G Z'^. If ^ denotes the 
lexicographic order in Z"^, define 

r = r{K,f3,h) = Ej^.yo (^E^^,j,,^olog(e^v^(^o-^«)-)^ J , 

where Jq is an independent copy of Jq, Jq being independent of {Ji,i G Z"^). 
Here, the E-expectations are conditional, and the subscripts of E indicate on 
which variables the expectation is performed. Then, we have the following 
central limit theorem: 
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Theorem 3. For < k < Ko{h) and < (3 < po{h), 

V^(pjv(«:,A/i; J)-Ep^(«:,/3,/i; J)) AT (^0, T - ^g^^ . (6) 

One can check by expansion around k = f3 = 0, h ^ 0, that the hmiting 
variance T — ^<f is strictly positive for any fixed non-zero h and small k, (3. 
This proves that the fiuctuations of pN{i<-,P,h] J) are truly of order of the 
square root of the volume inverse in this region of the parameters. To match 
with the breakdown in the order of magnitude of fiuctuations at zero external 
field, we observe that both q and F vanish as — > 0, and also that they are 
equal to zero when h = 0. 

3 Proofs 

3.1 For small k, the RFIM is inside the Dobrushin uniqueness 
region 

Let i khe lattice points. Under any (infinite volume) RFIM Gibbs measure 
(■)oo J oo J '^^ '^i given aj = rij,j ^ i, is Bernoulli with parameter 

proportional to exp{cri[7ij ^..^j + K,K{i — k)rik]} with 

Following e.g. Section 2 in we define the Dobrushin's infiuence coefficient 

Cfc, = sup I {■\v)i,oc,J - (■W)i,oc,J \Ur;V = V' off fcj 

where || ■ ||var is the variation norm. With a straightforward computation, 

II (■|^)*,oo,J- (■|^')i,oo,jllvar = i,oo,J ' i+W) i,oo,J \ 

^ 2\smh2KK{i - k)\ 

~ cosh[2?ij^fc;^] + cosh.[2KK{i — k)] 

for such 7], r]', and so 

Cki < I tanh 2KK{i - k)\ < \2KK{i - k) \ (7) 

Therefore, for k < Ki = (2 X]j^o 1-^(01)^^ derive from and ^ that 
a = supj ^j^. pl*~^ICjfcj < 1 for some p > 1, which implies (see Section 2.3 in 
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that the Gibbs measure i-)^ j is unique and has exponentially decreasing 
correlations. More precisely, for all local functions /, g, there is a finite 
constant C = C{a, f, g) such that for all i G Z'^ 

|(/;^7oft,)^_,|<Cp~l^ (8) 

with 9i the shift of vector i and (/; g)^ j the covariance of /, g. 

3.2 Uniqueness of q 
Introduce 

n) = Fp,,{q, k) = h, P^) + ^(1 - qf 

which is, in view of (jHl), a smooth function of all its arguments if k < ki (e.g.. 
Cor. 8.37 in P). 

Proposition 1. Assume h 0. For all j3, there exists k,2 = i^2{l3,h) > 
such that q \—>- F{q, k) has a unique minimizer on [0, 1] for n < k,2- 

By [7j, we know that the minimizer go of F{q, 0) is unique, strictly positive 
and that^ 

0^/(^0,0) >0. (9) 

By (jHl), the function F is continuous in k uniformly in g G (0, 1]. Hence, 

V(5 > 3k : \/k < k argmin^ F(-, k) C (go — 5,10 + ^) (10) 

Again by (jH)), the function F is in a neighborhood of (go,0). By the 
implicit function theorem for the equation 

(which applies thanks to condition (jHl) ) we then derive that there exist neigh- 
borhoods U of go and V of k = and a function q : V ^ U such that, for 
q E V, K E U, the above equation is equivalent to g = qi^n). With 6 small 
enough so that (go — 5, go + 5) C V, we choose now K2 < k (with k from (fTIH)) 
such that (— /«2, K2) C U. Then, for k, < H2, F{-,k,) has a unique minimum 
at g = q{n). 

As for the case h = 0, one can prove similarly the following results, that 
we mention for comparison but will neither prove nor use. 

p.166 
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Proposition 2. For all (3^1, there exists = i^siP) > such that q ^ 
-^/3,o(?5 k) has a unique minimizer on [0, 1] for n < k^. If P < 1, the minimizer 
is q = 0. 

Remark The restriction /3 7^ 1 is due to the fact that for /3 = /3c = 1 (the 
critical point of the Sherrington-Kirkpatrick model) one has qo{l3c,0) = 
and, in contrast with 

^F,,,o(0,0) = 0. 
3.3 Proof of Theorem Q] 

Proof of ^ The proof is standard, we just sketch the main steps. Consider 
the box AmN, with m,N G N, and partition it into sub-boxes A^'*, i = 
1, • • • ,m'^, congruent to A^. Moreover, let Z^n be the partition function of 
the system where, with respect to (H)), H^^{a; J) is replaced by 

the being m independent families of standard Gaussian variables. Note 
that, in this system, the different sub-boxes interact only through the Ising 
potential K{.). Then, following the ideas of JU]; it is easy to prove that 

PmN{i^,f3,h) > E log Z^n{i^, (3, h; /^)). 

Since the Ising potential is summable, the interaction among the different 
sub-boxes due to the potential K grows at most proportionally to k and to 
the the total surface, dm'^N'^^^ . As a consequence, one has the approximate 
monotonicity 

C 

PmN{K, 13, h) > pn{k, 13, h) - K — , 

for some constant C depending on the potential K{.). From this, it is a 
standard fact to deduce that the sequence {pmN}^ has a limit when N ^ 00, 
that it does not depend on m, and that it coincides with lim^vpTv 

Proof of ^ The almost sure convergence is standard and follows from 
exponential self- averaging of the pressure (see for instance Proposition 2.18 
of ^2]), which in the present case reads 

P {\pN{^^,P,h;J)-EpN{K,f3,h;J)\ > u) < g-^^^^)!'^-!"', 

together with Borel-Cantelli's lemma. The L^-convergence comes from uni- 
form integrability, which again follows from exponential concentration. 
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3.4 Proof of Theorem 

For < t < 1 and any q > 0, define the interpolating partition function 

Z{t)= Yl exp(-//W(a)) (11) 

cr(^{-l, + l}^N 

witli tfie properties 

ZiO) = Zff'^iK,h,P^; J) 
Z{1) = ZN{^i,P,h;J). 

Tlie t-derivative of the corresponding pressure 

PNit) = r^E log Z{t) 

I A AT I 

is easily computed: We denote by {■)^ the Gibbs measure associated to 
by its tensor product acting on a pair {a^,a'^) e {— Ij+l}'^^ x 
{ — Ij+l}^'^, by qi2 = \-^n\~^ '^i(zAj^ o'iO'f the overlap between configura- 
tions and a^, and we get by the Gaussian integration by parts formula 
EJiF{Ji)=EF'{Ji), 



int. by parts /? J (3 



^(i-(.>j.?4>f)-,3,5;(i-(.>?>f)l 

jeAjv isAjv J 



2|A^| 2|A^, 

2 q2 



= Y(l-#-YE((g,.-g)^>f (12) 

so that, integrating in t between and 1, taking the N ^ oo limit and 
optimizing on q, we have the "first half" of Eq. Q: 

p{K,P,h) < inf fp™(K,/i,/3v^) + ^(1 - qA . (13) 



q>Oy V ' - V 4 

Note that the inequality holds for any values of f3 and n. 

Next, consider a system of two coupled replicas: For A > let 

Z(^)(t, A) = Yl (-H^H^') - H^\<y') + y |A^|A(gi2 - qf) 
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and denote by Gibbs measure on { — x { — Ij+l}"^^ asso- 

ciated to this partition function. (Of course, (■)^q = and Z^'^\t,0) = 
Z{ty.) With 



2|A 



-E log Z^^\t,X) 



N\ 



we have by a computation similar to p2|l and using symmetry in cr^, o"^, 
I pg)(t,Ao-t) = 

-/?|A^|((gi2-g)^>,,„_; 



int. by parts 



/5 



/3 



4|Aiv| |A 



TV 



^ (1 + (<.,V>.V|>, - 2 (a.V>.V|>fI_,) 



-2/3,^(l + (..V.'>,,„_,-2K'< 

ieAiv 



,Xo~t 



^(l-g)^-yE((g,3-g)r' 



t,\o-t 



so that, integrating 

pg)(t,A)<^(l-gft + 



, 1 E log ^(^^(O, t + A). 



4 



2|A 



N\ 



(14) 



Now, setting MAr(t) = pff™{K,, h, l3y/q) + ^(1 — q)H — PN{t), which is a non 
negative function by (|T^ . and using convexity of the pressure with respect 

('2') 

to A and the identity Pj^ it, 0) = PN{t), we obtain for any A > 



d 1121 U (2). , 



d 

T) 

.(2) 



convexity {t, X) - p^it) 



< 

d 1 

< - 

A 



A 

UN{t) - 



(15) 



2|A 



1 Z(^){0,t + X) 

-^'""^ Z(2)(0,0) 



7V| 



Since Z*^^^ is increasing in A, Eq. ()15|1 implies 



2|A 



TV 



Z(2)(0,0) 



1 
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which, recalhng that un{0) = 0, can be immediately integrated to give 

and Eq. (0) follows if we can prove that 

lim -T^E log ^Zi^' = (17) 
iv^oo2|A^| ^ Z(2)(0,0) ^ ^ 

for some Aq > 1, if g is chosen properly. (Note that Aq = A + 1 > 1 is required 
so that one can take t up to 1 and still have A > 0, which is needed in (fT3j) .) 
Define for any /i G M, g > 

and a^ifJ') = ^C(N{f^] J)- We denote by jv Gibbs measure for two 

replicas of the RFIM with parameters k and 7 = (3^/q and volume A^r. Later 
we will also use the notation ^ for 



52 



/ /. M|Aiv|{gi2-g)\»- 

Let = gAr(/5, k, /;,) be the value which minimizes 

p™(«:,/^,/3y^) + ^(l-g)2 
with respect to g, cf Clearly, satisfies the "self-consistent equation" 

9 = E (912),,^^,^ = . (18) 

An analysis analogous to the one of Section 13.21 shows that the solution of 
(fTHj) is unique for k small enough, for sufficiently large, and that Qn Q 
for A^ — 00. A Taylor expansion around fi = gives immediately 

rt^ ry 

aN^fi; J) = aiy{0] J) + fia'^{0] J) + dy dua'l^{u] J), 
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where 

a^(0;J) = (19) 
«^(0;^) = 1{{Qu)Zv^,n-^{9i2)Zv^,n) (20) 

= ^ E ((-^-^-(-^-^»(-K-(-K»>!::Jv^,.(22) 

We can view (■)i"^^/^Ar as the Gibbs measure of a system with 2|AAr| spins, 
with exponentially decaying pair interactions. Since k is small, taking fi 
itself sufficiently small, we keep this system inside the Dobrushin uniqueness 
region. Using exponential decay of correlations as in Section 13.11 we obtain 
that a'l^{u] J) is bounded above by a constant, uniformly in N,J and in 
u G [0, fi], so that 

aAr(/i; J) < /ia^(0; J) + C/i^ (23) 

for K < Ko{h). Note that this bound holds for any /x, since it does for small 
II and the function a can grow at most linearly at infinity. 

With this in hand, we go back to proving (|T7jl . After a Gaussian trans- 
formation 



dz 



Z(2)(0,Ao) f , /|Arl .. / „ _ .X ®2 



Z(2)(0,0) J 



27r 



Equations ^ and {211 imply that, for (3 < (3o{h) = l/^/ACXo 



-^E log ^^ly ' < ^Elog/rf^J^el^-lHV2(i-4C/3^Ao)+2.,v/A^c 
2|A^| ^ Z(2)(0,0) - 2|A^| V 27r 

< C'E{a'^{0;J)y 

= J&r2 E E(((a,)^-E(a,)^)((a,)^-E(a/)), 

where for simplicity we have written (.) for (.)^2^^^^ \y"g show that 
the last expression is of order l/|Ajv|. Indeed, take two distinct sites i,j 
and consider the (i-dimensional ball B^j of radius \i — centered at i. If 
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(cTj) were depending only on {Jk,k G Bij) and (cXj) on {Jk,k G B^j) only, 
the corresponding term in would be zero by independence. But in the 
Dobrushin region, we can approximate (cTj) by the expectation of for the 
finite volume RFIM on Bij with an error which is exponentially small in the 
radius \i—j\/2, uniformly in the J^'s. Doing similarly with {(Jj), we conclude 
that 

for suitable constants C',C" > 0. This, together with (j21I), immediately 
implies that 

1 ^, z^Ho^^ c 



2\An\ " Z(2)(0,0) - lA^vl' 
At this point, recalhng Eq. (fTT)|l . one finds 

Pn{^, P, h) > pT^^k, h, PV^) + ^(1 - qr^f - (e^/^^^-^) - l) 



|A 



N\ 



which, together with (jlHj) . proves the convergence in average of the pressure 
to the expression (jSJ. Moreover, from Eqs. (fT3|) and (fT^ one deduces that 

limsup sup |Ajv|IE {{qu — qNT)f'^ < oo , (25) 

N^oo 0<t<l 

which will be needed in next section, where we deal with pressure fluctua- 
tions. 



3.5 Proof of Theoreml^ 

We follow the strategy which was introduced in [TT], adapted to the present 
case where short-range interactions are also present. Let 

V \^^n\ 

where Z{t),pjq{t) were defined in Eqs. (ITT|) - (fT^ and it is understood that 
q is taken to be qjq = qiy{P,K,h) as in Eq. (fTS|) . to be distinguished from 
q = lim q^. 

We will prove that 

lim Ee'"^^(*) = exp ( - — {T - —tf)] (26) 

N^oo \ 2 2 J 

for any mGM, 0<t<l, from which Eq. (jH)) follows for t = 1. 
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By means of integrations by parts one finds 

= ^u'qlEe^-M'^ - ^u'Ee^-Mt) ((^^^ _ ^^^^^^ 

and, using Eq. (j^SJ, 

^tEe^^^^W = ^l^^Ee^"^^(*) + o(l). 
Integrating in t, Eq. is then proven provided that we show that 

\im Ee^"^''^^") = exp (-yr^ , (27) 



i.e., a central hmit theorem for pressure fluctuations of the RFIM at high 
temperature. 

To this purpose, we employ the central limit theorem for martingales JS] , 
which we recall for convenience. Let {Q, A, P) be a probability space, 
{^k}o<k<n a filtration of A, for n e N, such that J^l^ = {0, fi}, and ^("^ = 
{^n,k}i<k<n a sequence of random variables adapted to JF^"). We denote 
by E"-''' (respectively P"'^) the expectation (respectively the probability) 
conditioned to JF^ and by V^'^ the conditional variance 

of a random variable X. We say that the triangular array {'Cn,fc}^>o i<fc<n 
asymptotically negligible if for any £ > 



5^P"''=-Mlen,.|>^)^0 (28) 



fc=i 

when n — s> oo. Then, the following holds: 



Theorem 4. Lei {^n,fc}„>o i<fe<n asymptotically negligible triangular 

array of square integrable random variables, and assume that for some F > 0, 

n 

^i?-'^-i(e.,fc)^0 (29) 
fe=i 
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and 

n 

E^"''"'(^->^)^r, (30) 



fe=i 

for n CO. Then, 



law 
k=l 

To simplify notations in our case, let |AAr| = n, h + P^/qNJi = hi and 
Pn,/, = P™(«:, h,(3,/^- J). 

Introducing the usual lexicographic ordering of the sites in A^v, we define 
J-"^, for A; G Ajy, as the a-algebra generated by the random fields Jj for 
i G An, i ^ k, and 

SO that 

,/\I^\ip™{K,h,PV^-J)-Ep§^'''{K,h,PV^-J))= J2^r.,k. (31) 

Of course, E'^'^{.) = Ej^^£gAjv/;-A:(-)5 by convention JF^ = {0, f2} if k 
precedes the first site in A^v- One can rewrite 



Jn ^ \ / n,h 



where is an independent copy of Jk - independent of (J;, / G Z'') -, and 
h'^ = h + l3^/q^Jl, so that 

\Uk\<n-^'^^.r^\K-hk\<Cn^'Whk\+C) 

for some constant C and 

This implies asymptotic negligibility (j^Hj) . since 

^'""^"^ (i4n,fc rj < - 
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In order to apply Theorem |3J we have to check conditions (j^^ and ()30|) . The 
first one is evident, since 



identically. As for the second, notice that 



V 



n,k—l 



-Ej,. I Ej, 



n 



c log ( e 



n,h 



(32) 



Let k correspond to a site "in the bulk" of Aat, i.e., assume that the distance 
between k and the boundary of A^r is larger than, say, r2^/*^^^\ In this case, 
we will write k G B^. We want to replace {■)nh with the unique 

infinite-volume Gibbs measure i-)^!^- To this purpose, note preliminarily 
that 



log 



n,h 



co,h 



< 2\h'^ — hk\. 



Moreover, thanks to Dobrushin's theorem, 



\ /n,h 

/g(?i^-?Jfc)o-fc\ 

\ / oo,h 



1 H m7 — 1. N_ \ Sn,h 



I oo,h 



with some (5„ h such that 



(33) 



lim sup sup|5„/i| = lim 



Denoting by An^k the event 



A, 



n,k 



\hi - hk\ < - log ( — 



and using the fact that | log(l + a;)| < D\x\ for |x| < 1/2 for some finite 
constant D, one can write 



log- 



In.h 



\ I oo,h 

Therefore one has 

1 



n,k 



V 



n,fc— 1 



n 

1 

n 



Ej, ( E.,, 



kh) + oil)) 



log(e 



oo,h 



+ -o 1 

n 
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where o(l) —>■ for n oo uniformly in /c G B^, 9^ is the shift of vector k 
and 



the subscript o referring of course to the origin of the lattice l/". Note that 
there is a residual n-dependence in 0, since the fields hi are defined through 
gAT, but this dependence is easily seen to be harmless thanks to the exponen- 
tial decay of correlations inside Dobrushin's uniqueness region, and to the 
fact that qN ~^ <!■ Finally, defining hi = h + P^/qJ^, the ergodic theorem 
implies that, for almost every J, 



where we used the fact that the contribution of the spins k ^ Bj^ vanishes 
for lA^rl — >■ cxo. At this point, all the conditions necessary to apply Theorem 
lUare fulfilled and, recalling (jHlll . one has ipTj) . which concludes the proof of 
Theorem IHl 
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